Triangular flow in heavy ion collisions in a multiphase transport model 
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We have obtained a new set of parameters in a multiphase transport (AMPT) model that are 
able to describe both the charged particle multiplicity density and elliptic flow measured in Au+Au 
collisions at center of mass energy y^sjvjv = 200 GeV at the Relativistic Heavy Ion CoUider (RHIC), 
although they still give somewhat softer transverse momentum spectra. We then use the model to 
predict the triangular flow due to fluctuations in the initial collision geometry and study its effect 
relative to those from other harmonic components of anisotropic flows on the di-hadron azimuthal 
correlations in both central and mid-central collisions. 

PACS numbers: 25.75.-q, 12.38.Mh, 24.10.Lx 



I. INTRODUCTION 

Studies of anisotropic azimuthal flows in heavy ion 
collisions at RHIC have provided important informa- 
tion on the properties of produced quark-gluon plasma 
(QGP) [J-Q- In particular, the large elliptic flow ob- 
served in experiments has led to the conclusion that the 
produced quark-gluon plasma is strongly interacting as it 
can only be explained in the hydrodynamic model with a 
very small viscosity or in the transport model with 

parton scattering cross sections much larger than those 
given by the perturbative QCD d [T^. With the large 
parton scattering cross section, the transport model has 
also been able to describe the hexadecupole flow mea- 
sured at RHIC More recently, the importance of 
the triangular flow, which originates from fluctuations 
in the initial collision geometry 
out in Ref. Unlike the elliptic flow. 



has been pointed 
the triangu- 
lar flow is less sensitive to the centrality or the impact 
parameter of the collision (l3l - [T5| . Also, a study based 
on the viscous hydrodynamics has shown that the vis- 
cosity in the quark-gluon plasma has a larger effect on 
the triangular flow than the elliptic flow in relativistic 
heavy ion collisions [ig} . Furthermore, it was suggested 
in Ref. [l^ and later shown in Ref. in a multiphase 
transport (AMPT) model, which includes both initial 
partonic and final hadronic scatterings, that the trian- 
gular flow may play an important role in the away-side 
double-peak structure seen in the di-hadron azimuthal 
correlations at RHIC. In the present paper, we extend 
the study of Ref. to investigate more quantitatively 
the triangular flow in heavy ion collisions at RHIC by ad- 
justing the parameters in the AMPT model, particularly 
the parton scattering cross section, to fit more recent 
experimental data such as the charged particle multiplic- 
ity density and transverse momentum spectra as well as 
their elliptic flow. We find that the resulting tri ang ular 
flow has smaller values than those shown in Ref. [I7| but 
still has an appreciable effect on the di-hadron azimuthal 
correlations as found in Ref. [l7| . 

This paper is organized as follows. In Sec. [Ill we re- 



view the AMPT model and discuss the parameters in 
the model that are relevant to the present study. In 
Sec. mil we describe the results on the charged parti- 
cle multiplicity density and transverse momentum spec- 
tra in heavy ion collisions at RHIC and their compar- 
isons with experimental data. We then show in Sec. IIVI 
the calculated charged particle elliptic flow in compari- 
son with the experimental measurements as well as the 
predicted triangular flow. In Sec. El we study the effect 
of anisotropic flows, particularly that of triangular flow, 
on the di-hadron azimuthal correlations. In Sec. I VII we 
discuss the specific viscosity in the initial partonic mat- 
ter produced in heavy ion collisions. Finally, a summary 
is given in Sec. IVIIl 



II. THE AMPT MODEL 

The AMPT model is a hybrid model [l^ with the ini- 
tial particle distributions generated by the heavy ion jet 
interaction generator (HIJING) model [l^- In the de- 
fault version, the jet quenching in the HIJING model 
is replaced in the AMPT model by explicitly taking 
into account the scattering of mini-jet partons via the 
Zhang's parton cascade (ZPC) model [2^. These par- 
tons are recombined with their parent strings after the 
scattering, which are then converted to hadrons using 
the Lund string fragmentation model. In the version of 
string melting, all hadrons produced from the string frag- 
mentation in the HIJING model are converted to their 
valence quarks and antiquarks, whose evolution in time 
and space is modeled by the ZPC model. After the end of 
their scatterings, quarks and aniquarks are converted to 
hadrons via a spatial coalescence model. In both versions 
of the AMPT model, the scatterings among hadrons are 
described by a relativistic transport (ART) model [2l| . 

In previous studies, it was found that the multiplic- 
ity of charged particles measured in heavy ion colli- 
sions at RHIC could be well described by the default 
AMPT model with modified values of a and b [12] in 
the Lund string fragmentation function f{z) cx z~^{l — 
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TABLE I: Values of parameters for the Lund string frag- 
mentation and the parton scattering cross section in previous 
studies (A) and in the present work (B). 



a 


b (GeV-^) 






A 2.2 


0.5 


0.47 


1.8 


B 0.5 


0.9 


0.33 


3.2 



z)"" exp{—b m\/z), where z is the hght-cone momen- 
tum fraction of the produced hadron of transverse mass 
mx with respect to that of the fragmenting string. 
To describe the measured elhptic flow required, on the 
other hand, the AMPT model with string melting to- 
gether with a larger parton scattering cross section a k, 
97ra^/(2yLt^), where as is the QCD coupling constant and 
/Lt is the screening mass of a gluon in the QGP, than that 
given in the perturbative QCD. Values for these param- 
eters ar e gi ven in the first row (A) of Table HI In these 
studies [U , the elliptic flow was calculated with re- 
spect to the theoretical reaction plane, corresponding to 
the second-order event plane ^2 = 0. As a result, this 
might have underestimated the elliptic flow that is deter- 
mined with non-zero ^'2 to take into account the fluctu- 
ations in the initial collision geometry as in the experi- 
mental analysis. We have recently shown, however, that 
both the charged particle multiplicity density and ellip- 
tic flow measured in heavy ion collisions at the Large 
Hadron Colhdcr (LHC) can be described by the AMPT 
model with string melting by using the default a and b 
parameters in the HIJING model together with a smaller 
QCD coupling constant and a larger screening mass as 
given in the second row (B) in Table U [IJ]. We note 
that the parton scattering cross sections for parameter 
sets A and B are about 10 mb and 1.5 mb, respectively. 
Although the latter has a smaller value, it has a more 
isotropic angular distribution. 



III. CHARGED PARTICLE MULTIPLICITY 
DENSITY AND TRANSVERSE MOMENTUM 
SPECTRA AT RHIC 

We first show in Fig. [T] the centrality dependence of 
the multiplicity density of mid-pseudorapidity charged 
particles in Au-|-Au collisions at ^Jsnn = 200 GeV ob- 
tained from the AMPT model with string melting for 
the two parameter sets in Table U and compare them 
with the experimental data from the BRAHMS Collabo- 
ration [2^ . For the relation between the centrality c and 
the impact parameter b, we use the empirical formula 
c = Trb^/tTin [2^ with the nucleus- nucleus total inelastic 
cross section a-^ ~ 705 fm^ calculated from the Glauber 
model. It is seen that the multiplicity densities from the 
parameter set A shown by filled triangles are larger than 
the BRAHMS data shown by filled squares at all central- 
ity bins, while those from the parameter set B, which are 
given by filled circles, are consistent with the BRAHMS 
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FIG. 1: (Color online) Centrality dependence of the multiplic- 
ity density of mid-pseudorapidity charged particles produced 
in Au-I-Au coUisions at = 200 GeV from the AMPT 

model with string melting using parameter sets A (filled tri- 
angles) and B (filled circles). The BRAHMS data shown by 
filled squares are from Ref. [25| . 
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FIG. 2: (Color online) Transverse momentum spectra of mid- 
pseudorapidity (|ry| < 0.18) charged particles by using param- 
eter sets A (filled triangles) and B (filled circles) in Au+Au 
coUisions at ./s^ = 200 GeV. The PHENIX data shown by 
filled squares are taken from Ref. [St! ]. 

Figure [2] shows the transverse momentum {pt) spectra 
of mid-pseudorapidity charged particles in different cen- 
trality bins (0 — 60%) from parameter sets A (filled trian- 
gles) and B (filled circles) . It is seen that both parameter 
sets describe reasonably the experimental data from the 
PHENIX Collaboration ^ at low pr- At high pr, the 
parameter set B gives a larger yield than the parameter 
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set A as a result of smaller energy loss of high-pT particles 
when the particle density is lower and the parton scatter- 
ing cross section is smaller as for the parameter set B. As 
in the previous study [l^ , the pr spectra from both pa- 
rameter sets are softer than the experimental data, and 
this is due to the small current quark masses used in the 
AMPT model so that partons are less affected by the 
radial flow effect. 



IV. CHARGED PARTICLE ANISOTROPIC 
FLOWS AT RHIC 



that for all centralities, the elliptic flow is larger for the 
parameter set A than for the parameter set B. This orig- 
inates from two effects. First, the larger string tension 
in parameter set A leads to a larger number of initial 
particles, which results in a larger pressure and thus a 
larger anisotropic flow. Second and more importantly, 
the larger parton scattering cross section in parameter set 
A converts more efficiently the initial spatial anisotropy 
to the final momentum anisotropy. As for the charged 
particle multiplicity density and transverse momentum 
spectra, the parameter set B gives a much better descrip- 
tion of the measured elliptic flow than the parameter set 
A. 




FIG. 3: (Color online) Transverse momentum dependence of 
the elliptic flow of mid-pseudorapidity {\rj\ < 0.35) charged 
particles by using parameter sets A (filled triangles) and B 
(filled circles) in Au-|-Au collisions at y^ijviv ~ 200 GeV from 
the two-particle cumulant method. The PHENIX data shown 
by filled squares are taken from Ref. [2^ . 

In Fig. [31 we compare the transverse momentum 
dependence of the elliptic flow of mid-pseudorapidity 
charged particles in different centrality bins (0 — 60%) 
of same collisions from the two parameter sets with the 
experimental data from the PHENIX Collaboration [2^ . 
In both theoretical calculations and experimental analy- 
ses, the elliptic flow is determined using the two-particle 
cumulant method [29l.[30j 

V2{2} = V(cos(2A</))), (1) 

where Acj) is the azimuthal angular difference between 
particle pairs within the same event and (• • • ) denotes 
average over all possible pairs. As the nonflow effect is 
overestimated at higher px in the AMPT model as shown 
in Ref. [131 1 '^^ only compare the anisotropic flows at 
lower pt where the nonflow effect is small. It is seen 



FIG. 4: (Color online) Transverse momentum dependence of 
the triangular flow of mid-pseudorapidity (I77I < 0.35) charged 
particles by using parameter sets A (filled triangles) and B 
(filled circles) in Au+Au collisions at ^/snn ~ 200 GeV from 
the two-particle cumulant method. 

We have also studied the transverse momentum de- 
pendence of the triangular flow of mid-pseudorapidity 
charged particles in different centrality bins (0 — 60%) 
in Au-I-Au colhsions at ^/snn — 200 GeV. Results 
obtained from the two-particle cumulant method, i.e., 
?;3{2} = a/ (cos(3A(/))), are shown in Fig. |4] using both 
parameter sets A (filled triangles) and B (filled circles). 
Compared with the elliptic flow shown in Fig. [3l the tri- 
angular flow is smaller and less dependent on the central- 
ity. Similar to the case of elliptic flow, the triangular flow 
is larger for the parameter set A than for the parameter 
set B. Since the parameter set B has been shown to give 
a better description of the charged particle multiplicity 
density, transverse momentum spectra, and elliptic flow, 
we believe that it would also give a more reliable pre- 
diction for the triangular flow. We note that the magni- 
tude of the triangular flow from the parameter set B is 
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similar to that from the (3 + l)d viscous hydrodynamic 
model [3 with the specific viscosity of the value 0.08 and 
the transport+hydrodynamics hybrid approach [isj but 
a little larger than that in Ref. ^ based on the AMPT 
model. 
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FIG. 6: (Color online) Comparison of the exact azimuthal 
correlations per trigger particle with the approximate ones 
defined in Eq. ([2]) up to order n = 5 in Au+Au collisions at 
^/sNN ~ 200 GeV for impact parameters b = fm ((a) and 
(b)) and b = 8 fm ((c) and (d)) from the parameter set B. 



FIG. 5: (Color online) Centrality dependence of the elliptic 
flow (filled symbols) and the triangular flow (open symbols) 
from parameter sets A (triangles) and B (circles) in Au+Au 
collisions at ^snn = 200 GeV for mid-pseudorapidity (|??| < 
0.35) charged particles in the transverse momentum window 
0.2 GeV/c < pt < 8 GeV/c from the two-particle cumulant 
method. 



Figure [5] displays the centrality dependence of the el- 
liptic flow (filled symbols) and the triangular flow (open 
symbols) of mid-pseudorapidity charged particles in the 
transverse momentum window 0.2 GeV/c < px < 8.0 
GeV/c from the two-particle cumulant method. Com- 
pared to the elliptic flow, the triangular flow shows a 
much weaker centrality dependence. As for the pt- 
dependent differential flows, the momentum-integrated 
flows are larger for the parameter set A than for the pa- 
rameter set B. Both elliptic and triangular flows increase 
with increasing centrality at small centralities but satu- 
rate or decrease at large centralities. Again, results from 
the parameter set B are similar to those from the trans- 
port+hydrodynamics hybrid approach 15 1. 



following approximate azimuthal correlations: 



DI-HADRON AZIMUTHAL CORRELATIONS 



The effect of anisotropic flows on the di-hadron az- 
imuthal correlations is a topic of great current inter- 
est [H, [13, IsTI - fssj and can be studied by considering 



dAcj) 



1 

+C50 



2 ^(ATt^SiVassoc^t^rig^assoc^^ COs(nA(^)], 

(2) 



n=l 



where (• • • )o denotes average over all events, N^"s and 
^assoc g^j.g numbers of trigger and associated particles, 
and and 1)^^°^^ are, respectively, their nth-order 

anisotropic flows. The nth-order anisotropic flow is cal- 
culated with respect to the corresponding event plane 
defined by 



1 (p^ sin(n(t>)) 

ipn = - arctan — ^ — — , 

n [p^ cos(n0)) 



(3) 



where (j) is the azimuthal angle and p^ is the weight fac- 
tor. As in Ref. [3l|, we take the transverse momenta of 
trigger particles to be p^'^ > 2.5 GeV/c, corresponding 
mostly from the fragmentation of energetic jets produced 
in initial hard scatterings that have not interacted much 
with the produced medium. For the associated parti- 
cles, their transverse momenta arc taken in the window 
1 GeV/c < p^''°'= < 2 GeV/c as in Ref. [H, and they are 
medium particles after the passage of initially produced 
back-to-back jet pairs. At mid-pseudorapidity {\r]\ < 1), 
the numbers of trigger particles and associated particles 
from the AMPT model are, respectively, 1.336 and 128.6 
for b = fm, and 0.565 and 42.3 for b = 8 fm. We first 
show in the upper panels of Fig. [H] by filled squares the 
exact di-hadron azimuthal correlations per trigger parti- 
cle calculated from all possible pairs of trigger particles 
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and associated particles from the AMPT model for the 
two impact parameters b = fm and b = 8 fm, cor- 
responding to central and mid-ccntral collisions, respec- 
tively. It is seen that for both impact parameters there is 
a peak around = at the near side of trigger parti- 
cles, while at their away side, i.e., around A(/) ~ tt, there 
is a broad structure for b = f m but a pronounced peak 
for b = 8 fm. As shown by dashed lines, the approximate 
di-hadron azimuthal correlations calculated with Eq. ([2]) 
using a = 1/4 in determining the event plane (Eq. ([3])) 
reproduce very well the exact ones, and their difference, 
shown in the lower panels of Fig. [51 is smaller than the 
residual correlations after subtracting the contributions 
from anisotropic flows as discussed below. 
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FIG. 7: (Color online) Di-hadron azimuthal correlations per 
trigger particle from anisotropic flows up to order n = 5 in 
Au-f Au collisions at ^/sNN ~ 200 GeV for impact parameters 
b = fm (left panel) and b = 8 fm (right panel) from the 
parameter set B. 



As in Rcf. [17|, we evaluate the contributions to di- 
hadron azimuthal correlations from anisotropic flows of 
various orders by replacing the event average of products 
with the product of event averages in the approximate 
correlations given in Eq. i.e.. 



dNr, 



dAcj) 

+ 00 



back 



-f2^(Ar'"%*"e)e(A^''""°'=wr°'>cCos(nA0)], (4) 

n=l 

where the first term is a constant independent of the 
azimuthal angle 0, and the other terms arc the contribu- 
tions from anisotropic fiows. In Fig. [71 we show the di- 
hadron azimuthal correlations per trigger particle from 
anisotropic flows up to order n = 5. It is seen that the 
elliptic fiow, which has a peak at A0 = tt, has the largest 
contribution for b = 8 fm, while for b = fm the con- 
tribution from the triangular fiow, which has peaks at 



Acf) = 27r/3 and 47r/3 in the di-hadron azimuthal correla- 
tions, is equally important. Furthermore, the triangular 
fiow is seen to give a larger contribution to the di-hadron 
azimuthal correlations per trigger particle for b = fm 
than for b = 8 fm as a result of the larger number of 
associated particles, although it has a larger value for 
b = 8 fm than for b fm as shown in Sec. IIVI As to 
the contributions from higher-order flows Vi and , they 
are, on the other hand, relatively small for both b = 
fm and b = 8 fm. 
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FIG. 8: (Color online) Di-hadron correlations per trigger par- 
ticle after subtracting background correlations up to different 
orders in Au-|-Au collisions at ^snn = 200 GeV for b = 
fm ((a), (b), (c) and (d)) and b = 8 fm ((e), (f), (g) and (h)) 
from the parameter set B. 



The residual di-hadron correlations per trigger particle 
after subtracting both the constant term in Eq. dH) and 
the contributions due to anisotropic fiows (Fig. [7]) from 
the exact azimuthal correlations (upper panels in Fig. [5]) 
are shown in Fig.[8l It is seen from Figs.[8ja) and (e) that 
the away-side double-peak structure is much stronger for 
b = fm than for b = 8 fm after subtracting the contri- 
bution from the direct flow as well as the large contribu- 
tion from the elliptic fiow. Subtracting also the contribu- 
tion from the triangular flow as shown in Figs, and 
(f) changes the away-side double peaks in the di-hadron 
azimuthal correlations to essentially fiat correlations for 
b = fm but to a single peak for b = 8 fm, while the 
near-side peak is reduced in both cases. As expected from 
the small contributions due to higher-order fiows shown 
in Fig. [71 further subtraction of the contributions from 
higher-order anisotropic fiows does not change much the 
shape of the di-hadron correlations as shown in Figs.[5Jc) 
and (g) as well as in Figs. [HJd) and (h). The residual 
di-hadron azimuthal correlations per trigger particle af- 
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tcr the subtraction of the contributions from anisotropic 
flows are then the correlations induced by initially pro- 
duced back-to-back jet pairs. Our results indicate that 
the away-side jets essentially disappear in the produced 
medium in collisions at b = fm but are still visible 
in collisions at b = 8 fm, consistent with the observed 
larger jet quenching in central than in mid-central colli- 
sions. The away-side double-peak structure in Figs, ^a.) 
and (e) after subtracting the contribution from the ellip- 
tic flow is thus from the combined effects of the away-side 
jets and the higher-order anisotropic flows, particularly 
the triangular flow. The weaker away-side double-peak 
structure for b = 8 fm in Fig. |8] (e) than for b = fm in 
Fig. [8] (a) then reflects the effect due to the remnant of 
away-side jets and the relatively smaller contribution per 
trigger particle from the triangular flow shown in Fig. [T] 
It is interesting to compare present results based on 
the parameter set B with those in Refs. and [3l| 
based on the parameter set A for b = 8 fm and b = 
fm, respectively. For the case of b = 8 fm, the away-side 
double-peak structure is weaker for the parameter set B 
than for the parameter set A, and the reasons for this are 
similar to those that cause its weakening in comparison 
to the case of b = fm based on the same parameter 
set. For the case of b = fm, the conclusion that the 
away-side double-peak structure in the results from the 
parameter set B is mainly due to the triangular flow is 
different from that in Ref. [3l| based on the parameter set 
A, which seems to indicate that other effects such as jet 
deflections or Mach cone shock waves are also relevant. 



VI. SPECIFIC VISCOSITY 



RHIC and LHC, it is of interest to estimate the specific 
viscosity in the partonic matter formed in these colli- 
sions. In the kinetic theory, the shear viscosity is given 
by r]s = 4(p)/(15(Ttr) in terms of the parton mean mo- 
mentum (p) and the parton transport or viscosity cross 
section atr = / dtd(j/dt{l — cos^6'). where t is the stan- 
dard Mandelstam variable for four-momentum transfer 
and da/dt k, Q'Ka1/^{t — /.t^)^] is the differential cross 
section used in the AMPT model. By assuming that the 
partonic matter only consists of non-interacting massless 
up and down quarks as in the AMPT model, we have 
(p) = ST and the entropy density s = (e -I- P)/T = 
4e/(3r) = 96T3/7r2 with T being the temperature of the 
partonic matter. The specific viscosity, i.e., the ratio 
between the shear viscosity and the entropy density, is 
then m 



' 40a^(9-f^)ln(i^i^)-18' 

The temperature dependence of the specific viscosity ob- 
tained from the parameter set B is shown in Fig. |9l and it 
shows that the specific viscosity decreases with increasing 
temperature. With the energy density of the baryon-free 
quark and antiquark matter given by e = 72T^ /tt"^, the 
initial temperature T is found to be about 378 McV in 
Au-I-Au collisions at ^/sWn = 200 GcV and 468 McV in 
Pb-|-Pb collisions at y/s^N = 2.76 TeV from the average 
energy density of mid-rapidity partons at their average 
formation time [l^. The specific viscosity is then about 
0.377 at RHIC energy and 0.273 at LHC energy as indi- 
cated, respectively, by the solid and open circles in Fig. [HI 
The values are thus similar at RHIC and LHC, although 
both are much larger than the lower bound p.^ / s ss 0.08 
predicted by the AdS/CFT correspondence [SJ]. 
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FIG. 9: (Color online) Temperature dependence of the specific 
viscosity in the partonic matter from the parameter set B. 

With the success of the parameter set B in describ- 
ing measured elliptic flow in heavy ion collisions at both 



VII. SUMMARY 

Using the default values for the parameters in the 
Lund string fragmentation function and a smaller but 
more isotropic parton scattering cross section than pre- 
viously used in the AMPT model for heavy ion collisions 
at RHIC, we have obtained a good description of both the 
charged particle multiplicity density and the elliptic flow 
measured in Au-|-Au colhsions at y^s^N = 200 GeV, al- 
though the transverse momentum spectra are still softer 
than the experimental results. With these constrained 
parameters, the magnitude of the triangular flow in these 
collisions has been predicted. We have also studied the 
di-hadron azimuthal correlations triggered by energetic 
hadrons at both impact parameters of b = and 8 fm 
and found that the double-peak structure at the away 
side of triggered particles, which is seen after subtracting 
the background contributions due to the elliptic flow, is 
largely due to the triangular flow. However, the residual 
correlations shown in our study after the subtraction of 
the flow contribution might still contain the contribution 
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from flow fluctuations besides the nonflow contribution 
that we are interested in [3^. It will be of great inter- 
est to find a method that can disentangle the nonflow 
contribution from that due to flow fluctuations. 

We have also estimated the specific viscosity in the 
initial partonic matter and found that it is much larger 
than the lower bound predicted by the AdS / CFT corre- 
spondence and is thus different from the values extracted 
with the viscous hydrodynamic model. The different con- 
clusions from the hydrodynamic model and the AMPT 
model might come from the fact that a constant specific 
viscosity is used in the former model while a constant 
total cross section is used in the latter one. Including 
the temperature dependence of the local screening mass 
in the evaluation of the parton scattering cross section in 



the AMPT model as in Rcf. [s^ may help to better un- 
derstand the different results from the transport model 
and the hydrodynamic model. 
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